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Abstract 

The winding of a single polymer in thermal equilibrium around a repulsive 
cylindrical obstacle is perhaps the simplest example of statistical mechanics 
in a multiply connected geometry. As shown by S.F. Edwards, this problem 
is closely related to the quantum mechanics of a charged particle interacting 
with a Aharonov-Bohm flux. In another development, Pollock and Ceperley 
have shown that boson world lines in 2 + 1 dimensions with periodic bound- 
ary conditions, regarded as ring polymers on a torus, have a mean square 
winding number given by (W 2 ) = 2n s h 2 /mk B T, where m is the boson mass 
and n s is the superfluid number density. Here, we review the mapping of the 
statistical mechanics of polymers with constraints onto quantum mechanics, 
and show that there is an interesting generalization of the Pollock-Ceperley 
result to directed polymer melts interacting with a repulsive rod of radius 
a. When translated into boson language, the mean square winding num- 
ber around the rod for a system of size R perpendicular to the rod reads 
(W 2 ) = 2™mk B T l n (^/°)- This result is directly applicable to vortices in Type 
II superconductors in the presence of columnar defects. An external current 
passing through the rod couples directly to the winding number in this case. 



I. INTRODUCTION 



The study of the statistical mechanics of polymers in multiply connected geometries 
began many years ago with work by S.F. Edwards and by S. Prager and H.L. Frisch 
The simplest nontrivial geometry consists of a polymer interacting with a repulsive rod, and 
the corresponding path integrals can be analyzed via an elegant analogy with the physics of 
a quantum mechanical particle interacting with a solenoidal vector potential || . The mean 
square winding number of the polymer around the rod, and other interesting quantities can 
be computed for this problem H. The physics bears a close mathematical relation to the 
famous Aharonov-Bohm effect for a real quantum mechanical particle interacting via its 
charge with a tube of magnetic flux ||. 

More recently, Pollock and Ceperley have studied the winding numbers with respect to 
periodic boundary conditions for boson world lines in the Feynman path integral formulation 
of superfluidity in two dimensions || . The physics here is equivalent to many ring polymers 
interacting on the surface of a torus. The mean square winding number of the world lines 
around the torus can be expressed exactly in terms of the renormalized superffuid density 
of the equivalent boson system 0. 

Winding numbers and the statistical mechanics of many polymer-like objects are also 
relevant to the physics of vortex lines in Type II superconductors ||. Here, thermal fluctu- 
ations in the trajectory of a vortex defect in the superconducting order parameter can be 
described by a Feynman path integral for an elastic string. A collection of many such lines 
behaves like a directed polymer melt, with the added complication of a quenched random 
disorder potential . Point-like disorder is always present to some degree due, for example, 
to quenched fluctuations in the concentration of oxygen vacancies in the high temperature 
cuprate superconductors. Drossel and Kardar have studied how the winding number distri- 
bution of a single directed polymer around various obstacles is affected by point-like disorder 

Dramatic improvements in the pinning efficiency of vortex lines have recently been 
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achieved via the introduction of columnar damage tracks created by heavy ion irradiation. 
If the concentration of damage tracks (assumed to pass completely through the sample) 
exceeds the number flux lines, there is a low-temperature "Bose glass" phase, in which every 



vortex is trapped on a columnar defect [[□]]. At high temperatures, the vortices delocalize 
in an entangled flux liquid. If these vortex trajectories are viewed as the world lines of 
quantum mechanical particles, the physics in thick samples becomes equivalent to that of a 
low temperature boson superfluid in 2 + 1 dimensions ||. 

In this paper, we discuss how to compute polymer winding numbers from the Hamil- 
tonian formulation of quantum mechanics, keeping in mind applications to vortex lines 
interacting with many columnar defects. A random distribution of columns maps onto a 



time-independent random potential in the quantum mechanical analogy We assume 
high temperatures, and samples which are clean before irradiation, so that point disorder 
is negligible and the usual Abrikosov flux lattice is either melted by thermal fluctuations or 
destroyed by the columnar disorder. An experimental realization of a multiply connected 
geometry for thermally excited flux lines is illustrated in Fig. 1. Superimposed on a dilute 
concentration of parallel columnar defects scattered randomly throughout a superconduct- 
ing sample is a slender tube in which the columns are very dense. Such a "tube of columns" 
could be made by covering a sample with a mask containing a small submicron hole during 
irradiation with a strong dose of heavy ions. Imagine that this sample is first subjected to 
a very large magnetic field, such that the density of flux lines is approximately equal to the 
density of columns in the tube. It should be possible to choose the temperature such that 
the vortices in the tube are in the Bose glass phase, while those outside constitute a flux 
liquid. The many trapped vortices inside the tube should then present a virtually impene- 
trable barrier to the thermally excited vortex lines outside. The concentration of lines in the 
liquid outside this repulsive cylindrical obstacle could be varied by decreasing the magnetic 
field. Because flux creep in the Bose glass phase is quite slow PJTl|, the concentration of 
vortices in the tube should remain approximately constant as a field is turned down, leaving 
the barrier almost unchanged. 



We shall explore the winding number fluctuations of the vortices with respect to the 
repulsive cylindrical tube in this experiment. As discussed below, a current passing through 
the tube couples directly to the net vortex winding number, and the mean square winding 
number of the unperturbed system gives the linear response of the net winding number to 
this current. (See Fig. 2) Such a current acts like an imaginary vector potential when this 
problem is mapped onto quantum mechanics jO]] . Winding of vortices around a thin repul- 



sive obstacle (or a set of such obstacles) could be probed via double-sided flux decorations 



12|| or indirectly by monitoring the flux flow resistivity in the plane perpendicular to the 
common direction of the applied field and the columns. The extra entanglements induced 
by the longitudinal current should impede vortex transport in this plane, and it would be 
especially interesting to look for changes in the in-plane resistivity as a function of the cur- 
rent through the tube [[HJ. This resistance should drop with increasing longitudinal current 



through the tubes, due to the enhanced vortex winding about the obstacle. 

In section II, we briefly review the Feynman path integral description of bosons in 2 + 1 
dimensions, and discuss the Pollack-Ceperley result for the winding numbers of bosons on a 
torus. We then describe the closely related statistical mechanics of directed polymer melts 
and vortex lines, and show how a current through a cylindrical obstacle couples to the 
winding number. In section III, we treat winding number statistics for isolated polymers 
both with and without columnar pin disorder. Results for many polymers or flux lines 
winding about an obstacle are presented in section IV. 

II. TWO-DIMENSIONAL BOSONS AND DIRECTED POLYMER MELTS 

A. Boson Statistical Mechanics 

The partition function for a set of N nonrelativistic bosons interacting with pair potential 
V(r) reads 

Z = Tr'ie-P*"), (2.1) 
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where (3 = 1/ksT and the boson Hamiltonian operator is 

N -h 2 l 

^ = E^rV? + -En|r,-r,|). (2.2) 

The prime on the trace means that only symmetrized boson eigenfunctions are to be included 
in the partition sum, and we shall focus on particles in two space dimensions. This trace can 
be rewritten in terms of a Feynman path integral by breaking up exp[— f37i] into M pieces 
(M»l), 



e 



-PHb 



e -ATH be -ArH b ... e -ArH b ^ ^3) 



where ArM = f3. Upon inserting complete sets of position states between various terms in 
the product and taking the limit M — > oo, the boson partition function may be expressed 
as an integral over a set of polymer-like trajectories {^-(r) = [xj(r), yj(r]} in imaginary time 
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1 N 
• P j=X 



2hp[Jo \dr) T 



1 rPK 



(2.4) 

/ii 

The normalized sum over permutations P insures that only boson eigenfunctions contribute 
to the sum, where the trajectories obey a type of periodic boundary condition, 

{r,-(/B»} = P[{r 4 (0)}], (2.5) 

and the operator P permutes the set of starting points {r,(0)}. 

An approximate picture of the polymer statistical mechanics problem represented by 
( p.4|) can be constructed as follows fi~4jl : Provided interaction effects are not too large, each 



'polymer" simply diffuses in the imaginary time variable r 

H 



r i (r)-r i (0)n«2-r (2.6) 



rn 



where h/m plays the role of a diffusion constant. When the temperature is high, only the 



identity permutation contributes to Eq. (|2.4j) . When projected onto the (x, y)-plane, the 



boson trajectories then behave like small ring polymers of typical transverse size given by 
the thermal deBroglie wavelength At, 

A T = 2nh 2 /mk B T 

~(|r,(r)-r,(0)| 2 )| T= ^. (2.7) 

At temperatures low enough so that At>u~ 1 / 2 , where n is the particle number density, 
complicated cyclic permutations appear, as the trajectories coalesce to form much larger 
rings. Feynman suggested in 1953 that the lambda transition from a normal bulk liquid of 
He 4 to a superfluid is associated with a dramatic proliferation in the number and length of 
such cooperative ring exchanges fL5| . 

Now, following Pollock and Ceperley, consider what happens when the excursions of 
the bosons in the xy-plane occur in a two-dimensional periodic box of size D |||7j (see 
Fig. 3). The permutation requirement ( |2.6[ ) describing periodic boundary conditions in the 
r-direction remains in effect. At high temperatures, virtually all ring polymers return to 
their initial positions {r-j(O)} when r = f3h, and the spatial periodic boundary conditions are 
unimportant. In the low temperature limit, however, ring exchanges lead to huge composite 
trajectories which typically wrap completely around the torus embodied in the (x, y)-plane 
periodic boundary conditions. The collection of cyclic boson trajectories can be classified 
by a dimensionless topological invariant, the vector winding number W, 

1 N 

W = ^=£M/%)-r,(0)] (2.8) 

V" j=i 

where Q = D 2 is the cross-sectional area of a square box with periodic spatial boundary con- 
ditions. In evaluating Eq. (|2.8|) , we imagine that the {^(t)} pass smoothly into neighboring 
periodic cells, without invoking the periodic boundary conditions (see Fig. 3). Consider the 
mean square winding number (W 2 ), where the angular brackets represent a path integral 
weighted by the exponential factor in Eq. ( |2.4| ) and divided by the boson partition func- 
tion. Ring polymers which do not wrap completely around the torus make no contribution 
to Eq. ( |2.8|) . However, at temperatures low enough so that >> n^ 1 ^ 2 , virtually all 
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trajectories belong to a cycle with a nontrivial winding number. We can then regard Eq. 
(p.8|) as a normalized iV-step random walk with typical step size Ay, and estimate that 



(W 2 ) w = 2-Knh 2 /rnksT, where n = N/Q is the number density of bosons. Pollock 



n 

and Ceperley showed that (W 2 ) is in fact exactly related to the superfluid density n s of the 
equivalent boson system, 

(W 2 ) = 2n s (T)h 2 /mk B T 

= 2n s (T)(h 2 /m)f3. (2.9) 

The implications of this remarkable connection between a topological invariant for a set 
of ring polymers on a torus and the superfluid density of the equivalent set of bosons is 
summarized in Fig. 4, where (W 2 ) is shown as a function of /3, which is proportional to 
the number of "monomers" if the trajectories are viewed as ring polymers. Translational 
invariance of the boson system in the absence of disorder implies that lirciT^o n s{T) = n, 
fl6fl , so we know that (W 2 ) diverges linearly with f3 as (3 — > oo, with slope given exactly by 
2nh 2 /m. The existence of a sharp Kosterlitz-Thouless phase transition in superfluid helium 
films, moreover, implies that there must be a singularity in (W 2 ): This quantity vanishes 
below a critical value f3 c , and the universal jump discontinuity in the superfluid density fL7| , 
lim T ^ T - n s (T)/T = I&bS, implies the exact result || 

lim (W 2 ) = -. (2.10) 

Feynman hoped that his "cyclic ring exchange" picture of superfluidity could be used to 
understand the lambda transition in He 4 which was a quite mysterious phenomenon in 1953 
|14] , [15| . Here, the sophisticated understanding of vortex unbinding transitions in He 4 films 
developed in the past 25 years (using other methods) has been used to make a prediction 
about a topological quantity, the winding number. 

Equation ( |2.9| ) allows winding number statistics extracted from computer simulations 
of bosons with periodic boundary conditions to be converted into measurements of the 
superfluid density . In can be used, in particular, to probe the reduction in the superfluid 



density of helium films due to substrate disorder near T = 0. A disordered substrate 
potential maps onto randomness correlated along the imaginary time direction r in the 
2+1-dimensional world line picture of boson physics. The winding number fluctuations are 
reduced because this correlated randomness reduces the wandering of the boson world lines. 



One might think that the results described above would be applicable, at least in princi- 
ple, to real ring polymers on the surface of a torus, whose fluctuations include fusing together 
to form large rings, similar to sulfer ring polymers in equilibrium fl8f . Unfortunately, the 
boson model as applied to ring polymers is unrealistic: The individual "monomers" within a 
ring, indexed by the imaginary time coordinate r in Eq. ( |2.4| ), are noninteracting, so there is 
no intra-chain self-avoidance. The only mterpolymer interactions, moreover, occur between 
monomers with same imaginary time coordinate, which is also unrealistic for a melt of ring 
polymers. 

The physics of directed polymer melts, on the other hand, is much closer to that of real 
bosons. Because the polymers are directed, interpolymer interactions at the same "imaginary 
time" coordinate dominate the physics. Vortex lines in superconductors provide an excellent 
physical realization of this system, but examples can also be found in polymer nematics in 
strong external magnetic or electric fields ]TjJ. Consider a collection of N vortex lines or 
polymers in three space dimensions, labelled by x, y and r (see Fig. 5). We assume that 
these lines are stretched out on average along the r axis, so that they can be described 
by single- valued trajectories {rj(r)}. The partition function is a multidimensional path 
integral, similar to Eq. Q2.4| ). 



B. Directed Polymer Melts 




3=1 



(2.11a) 



where 




(2.11b) 
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Here, the mass m in the boson problem has been replaced by a line tension g, % has been 
replaced by the temperature T and the polymer system thickness L plays the role of (3h 
for the bosons. Unless indicated otherwise, we shall henceforth use units such that ks — 1. 
V(r) is a repulsive interparticle potential and U(r) represents the random potential due 
to columnar pins (not shown in Fig. 5). Both V(r) and U(r) are independent of r. This 
approximation is quite reasonable for directed polymer melts and vortex lines in 2 + 1 



dimensions, provided (|;r| 2 ), the mean square tilt away from the z axis, is small fll 



The statistical mechanics of directed polymers differs from that for bosons due to the 
absence of periodic boundary conditions in the "imaginary time" direction r: we assume free 
boundary conditions in Eq. ( |2.11a| ), i.e., we integrate freely over the starting and end points 
of the polymers. Hence there is no sum over permutations and no condition analogous to 
Eq. (|2.5|) . This change, however, is less severe than one might think. Indeed, the directed 



polymer partition function can be written in a form similar to Eq. ( |2.1| ), || 

N . N , 

Z=Ii J mi dv j {v' l ...v' N \e-' HL ' T \v l ...v N ) (2.12) 

j=i i=i 

where 

rp2 N N 

n = E V? + E U(r 3 ) + E V(\ri - ril), (2.13) 

3=1 3=1 i>3 

and the states |iv--r/v) and (r^_ - - • r^y| describe entry and exit points for the polymers 
or vortices at the top and bottom of the system. The Hamiltonian 7i is identical to the 
Hamiltonian 7i b in Eq. ( |2.2|) , provided we introduce a disorder potential U(r) to model 
the effect of, say, a random substrate on the bosons. Because the initial and final states 
involve symmetric integrations over the entry and exit points, only boson eigenf unctions 
contribute to the statistical mechanics defined by Eq. ( |2.12| ). As, L — > oo, the physics will 
be dominated by a bosonic ground state and bosonic low- lying excitations, just as in Eq. 
( P-ip . Thus, the precise choice of boundary condition is irrelevant in the "thermodynamic 
limit" of large system sizes. 

Can one define meaningful "winding number" problems for directed polymers with free 
ends? Because we are no longer dealing with ring polymers, the "winding number" W in 



Eq. (|2.8|) (with f3h replaced by L) no longer has a precise topological meaning, even for 
periodic boundary conditions in the space-like directions. One can, however, still define and 
estimate the order of magnitude of the fluctuations in W as before. In terms of quantities 
appropriate for directed polymer melts, we have, up to constants of order unity 

(W 2 )^nTL/g, (2.14) 

where n is the areal number density of polymers in a constant r cross section. As illustrated 
(in the absence of a disorder potential) in Fig. 4, we expect that this quantity is nonzero 
for all values of L = (3% > 0. Although the curves for directed polymers and bosons should 
agree when L and f3h are large, there is no reason to expect a sharp Kosterlitz-Thouless 
phase transition for "bosons" with open boundary conditions |[20|| . The winding number 



variance for real bosons only vanishes identically for f3 < (3 C because the contribution to 
W from "ring polymer" trajectories which are not wrapped completely around the torus is 
exactly zero. There is no such condition for the open boundary conditions appropriate to a 
directed polymer melt. 

C. Winding Around an Obstacle 

A more interesting case is the winding number distribution for many directed polymers 
in a disorder potential and interacting with a repulsive cylindrical obstacle. As discussed in 
the Introduction, this problem has a direct physical realization in the context of thermally 
excited vortex lines. The combination of a disorder potential and repulsive barrier is illus- 
trated in Fig. 1. We model a repulsive obstacle by modifying the disorder potential in Eqs. 
( [HIED and ( gl5D , 



U(t)^U(t) + U b [t(t)\, (2.15) 



where the barrier potential Ub{t) (centered on the origin) is, e.g. 



oo, r < a 

U B (r) = { (2.16) 
0, otherwise. 
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This problem is a natural generalization to many directed lines of the original multiply 
connected geometry problem for single polymers of Edwards |I[ and of Prager and Frisch 
0. The (scalar) winding number with respect to the obstacle at the origin is now 

N r L dr- 
W = E/ ^A[r,(r)]-^, (2.17) 

where the "vector potential" A(r) is Q 

> , s t x r , 
A(r) = — 2 . (2.18) 

The boundary conditions along the r direction can be periodic or open ||. Periodic boundary 
conditions would imply integer winding numbers, but the noninteger winding numbers which 
characterize open boundary conditions are of direct physical interest for vortices, as we now 
explain. 

A current in the r direction through a cylindrical obstacle couples directly to the winding 
number of vortices outside in Type II superconductors. To see this, assume that the current 
is entirely confined to the cylinder, i.e., it cannot penetrate into the remainder of the sample. 
This "confined current approximation" should be applicable to a tube of columnar pins filled 
with pinned vortices in the bose glass phase, provided its conductivity is much higher than 
that of the flux liquid outside. Essentially infinite conductivity can be arranged in the tube 
because irradiation to produce columnar defects with matching field ~ 4 Tesla typically shifts 
the irreversibility temperature by 5-20°K. Choosing a temperature above the irreversibility 
line of the lightly irradiated region outside the tube but below the Bose glass transition inside 
insures that essentially all the current is confined to the tube. Under these conditions the 
additional free energy F of the vortices due the current is given by |21 
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6F = ~h f d2r J L(1tH ^ ■ B ^( r > r ) ( 2 - 19 ) 



where H ex t(r) is the field due to the extra current / through the wire 



, , 21 t x r . 
H cxt (r) = -— ^ (2.20) 



11 



and the prime on the integration means that the region inside the obstacle is excluded. The 
perpendicular component of the magnetic field Bi due to the vortices obeys the anisotropic 
London equation |2^ 



N rlv ■ M 

B ± (r, r) = O £ [r - r,(r)] + |^A 2 ViB ± (r, r) + A 2 «9 2 B ± (r, r), (2.21) 

j=i ± 

where 0o is the flux quantum, is the mass anisotropy and A is the screening length in the 
direction perpendicular to f. Upon substituting Eqs. Q2.20| ) and (|2.21| ) into (|2.19|) , integrat- 



ing by parts and noting that (H^A 2 Vj_ + A 2 <9 2 ) H ext (r) = in the domain of integration, 
we find that 

5F= JJoj : A^/xr, drj 



c f^[Jo 1i\r 2 j dr 
= -jW, (2.22) 

where we have defined a winding number coupling constant 7, 

7 = J0 o /c. (2.23) 

Consider now the mean winding number (W) induced by the current for vortices outside 
the obstacle. Upon incorporating an additional factor e~ SF ' T into the Boltzmann weight in 
( p.llap , we have, upon expanding (W) to leading order in the current, 

(W) = ^(W 2 )\ I=0 I. (2.24) 

Thus the equilibrium winding number fluctuations (TU 2 )|i=o determine the linear response 
of the net winding number to an external current. 

III. WINDING NUMBER AND QUANTUM MECHANICS: SINGLE POLYMER 

IN A CYLINDRICAL SHELL 

We now study the winding number fluctuations of a single directed polymer, moving on 
average along the r axis, and interacting with a repulsive cylindrical obstacle of radius a 
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centered on the origin of an (x, y) plane perpendicular to f. As discussed above, many of the 
results are directly applicable to vortex filaments in superconductors. We use the original 
method of Edwards which maps the generating function for winding numbers onto the 
quantum mechanics of particles interacting with a solenoidal vector potential. Although 
winding of a single polymer is simple and relatively well understood [P^JTm, we review the 
basic results to illustrate the method and discuss effects due to columnar pins outside the 
obstacle. We use the directed polymer notation discussed at the end of section II because the 
neglect of self-avoidance is justified in this case. The single line results would, however, apply 
in principle to a "phantom" polymer without self-avoidance in two dimensions interacting 
with a repulsive disk and a random substrate potential. We defer new results for many 
directed polymers to section IV. We shall, for simplicity, usually impose periodic boundary 
conditions in the r direction. This choice should not affect the results for (W 2 ) in the limit 
of large sample thicknesses L. For a systematic study of the effects of different boundary 
conditions in the context of the boson mapping, see ref. . 



A. Winding Number Formalism 

Consider the winding number ( p,17|) , specialized to the case of a single directed polymer 
with trajectory r(r), 

W= [ L drA[r(r)] ■ r(r) (3.1) 
Jo 

with A(r) = T^pz and r(r) = dr/dr. Note that A(r) is the electro-magnetic vector potential 
that would have been present had we replaced the cylindrical obstacle by a solenoid enclosing 
one unit of magnetic flux. If we then regard r as time, our expression for W becomes 
the quantum mechanical Aharonov-Bohm phase accumulated by a particle of unit charge 
traversing a path r(r), due to its interaction with the solenoid. As noted by Aharonov and 
Bohm |J, this phase is indeed the winding number, which is an integer for periodic boundary 
conditions, but can assume arbitrary values otherwise. 
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Our interest is, of course, not in the winding number corresponding to a single trajectory, 
but rather in the statistics of winding numbers. More concretely, we wish to study the 
distribution of W, averaged over all paths, with the appropriate Boltzmann weight. The 
winding number variance for a single polymer with line tension g is 



JVr(r) 



(W 2 ) 



where the free energy density 



JdT[i.A(v(r)] 2 e-Io^ir(r),r(r)] d r/T 



1 dr\ rT . . 



)}dr/T 



(3.2) 



(3.3) 



plays the role of a Lagrangian in the imaginary time path integral formulation of quantum 
mechanics. Here, U(r) includes the potential due to the columnar (i.e. r-independent) 
disorder and a part representing the excluded volume interaction with the cylinder, as in 
Eq. fl2JL5|). 

Upon defining a new function 



dr 

CMr),r(r)] = C -zTu-.A[r(r)}, 



(3.4) 



we have 



(W 2 ) 



d 2 



dv 2 



InZ(v) 



u=Q, 



(3.5) 



where 



T>t(t) exp 



1 r L 



TJo 



C u [r(r),r(T)]dT 



(3.6) 



It is straightforward using standard path integral manipulations to express Z{y) in terms 
of H(v), the Hamiltonian operator associated with C u , 



H{y) 



2g 



-V + uA(t) 



+ U(t) 



(3.7) 



specifically, 



Z(v)=Tr{e^WT} 



(3.8) 
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where the trace means we have now imposed periodic boundary conditions. For a system 
with cross-sectional area Q the mean square winding number itself is 

L d 



(W 2 ) = ~ J dviUr)) ■ A(r)U 



(W 2 ) = - 



L l d 2 H{v) 
T \ d 2 v 



v=0 



on, .an, , 
^7 (r) ^7 (r 



(3.9a) 



(3.9b) 



v=0 



where 



1 dH T 

9 



TdA 



-V + z/A(r) 



(3.10) 



is the current operator, (O) = itr {Oe T is the time ordering operator, and time 

dependence of operators is defined according to the Heisenberg picture in imaginary time 



0(i 



-Hr 



Oe Hl 



Equations ( |3.9a| , [3ljb| ) are central to our discussion. On the left-hand side they both 
have the fluctuations in the winding number. On the right-hand side, the first equation 
has the derivative of the persistent, thermodynamic, current due to a fictitious quantum 
particle flowing in a multiply connected geometry threaded by an Aharonov-Bohm flux. 
The derivative is taken with respect to the particle "charge" u, and the current is weighted 
by the vector potential. Equation ( |3.9b| ) expresses this derivative in terms of one- and two- 
time particle correlation functions. An analogous equation for many directed polymers in 
terms of boson correlation functions will be given in section IV. Information known from 
studies of quantum mechanical particles and many particle boson systems regarding these 
correlation functions will allow us to determine (W 2 ). 

Equations (|3.9a|J3.9b|) are expressions for the variance of the winding number W. We 
now briefly discuss the probability that the winding number W equals n, denoted by Pwin). 
This probability is, 

V w (n) 



Pwin) 



(3.11) 



where 
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Twin) = J Dr(z)5 J drA[r(z)] -r-n 



■ J Q L C [v(z),v(z)]dz/T 



r(r)]-r— inn 



2tt 



dvZ{y)e 



(3.12) 



with Z{y) given by (|3.6|). Equation ( |3.11 ) is again a mapping of the statistical property 
we are interested in onto a quantum mechanical problem, albeit a less transparent one: 
the probability Pw{ n ) is proportional to the Fourier transform, with respect to u, of the 
quantum mechanical partition function of a particle of charge v on an Aharonov-Bohm ring 
threaded by a flux 1/2tt. For a related discussion, including an earlier derivation of Eq. 
([nf), see Ref. [||. 



The normalization of the probability distribution function Vw{ n ) is a matter of some 
subtlety. In the second line of Eq. Q3.12|) we used the identity, 



drA[r(r)] ■ r — n 



2tt 



due 



drA[r(r)]-r— n 



(3.13) 



The limits of the v integration are determined by the allowed values of the winding number 
/ drA[r(r)] ■ r. Since we assume periodic boundary conditions, the allowed trajectories are 
closed (r(0) = r(L)), the winding number is an integer, and the integral over v can be taken 
between and 2n. However, if the allowed trajectories are not closed, the winding number 
can be non-integer, and the integral over v should be taken between — oo and +oo. For 
periodic boundary conditions, the integration domain v 6 [0, 2n] leads to J2 n 'Pw( n ) = 1- 
Although we have focused on periodic boundary conditions, in the limit L — > oo one expects 
the difference between periodic and open trajectories to vanish. 



B. A Single Polymer on a Thin Cylinder 

We now specialize to a polymer, whose motion in the x — y plane is confined to a ring of 
radius a, but is stretched out along the r axis. The polymer's Lagrangian is given by Eq. 
( [3 .3D and we add to U (r) a piece which confines the particle to an annulus just outside the 
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cylindrical barrier of radius a. Our previous discussion maps its random walk around the 
ring onto the thermodynamics of a quantum charged particle on an Aharonov-Bohm ring. 
The solenoid's vector potential is then A = ^— = constant, where is a unit vector around 
the ring. 

In the absence of disorder, we easily reproduce the expected results. According to Eq. 
( |3.9b|) the winding number variance is 

TT 

W) = —— 2 [l-2L(E{v = m (3.14) 

where (E(y = 0)) is the average kinetic energy of a chargeless particle on a ring. For a very 
large L the average energy is, to exponential accuracy, the energy of the ground state, i.e., 
0, and (W 2 ) = ( 2 nl) 2 g • ^ s ex P ec ted for a random walk, the variance of the winding number 
is proportional to the number of steps taken (i.e., proportional to L). The probability 



distribution P\y(n) is a Gaussian in this limit, as expected [10]. When L — > 0, (E) — > 
and (W 2 ) -> 0. 

The presence of a disorder potential U(r) modifies these results in an interesting way. 
The effect of disorder on 

^ = f \dl) 

—2iia d , „, . „ . 

lnZ(i/), (3.15) 



L dv 

i.e., the persistent current in a mesoscopic ring threaded by an Aharonov-Bohm flux, is 
well known: the disorder introduces a length scale, £, the quantum mechanical localization 
length of the ground state wave function. As L — > oo, 

Z(y)Kexp[-Leo(v)/T\, (3.16) 

where Eq(u) is the ground state energy. As a crude, but illuminating model of localization, 
consider a single narrow trap of depth Uq running up the side of the cylinder. The periodic 
boundary conditions around the cylinder then lead to a tight-binding model type result for 
the ground state energy as a function of u, 
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e (u) = -U Q - 2ti cos(u) - 2t 2 cos(2z/) (3.17) 

where the couplings t n represent the matrix elements for tunneling around the cylinder n 
times, t n ~ e - 27ran /C_ p or a square well of depth Uq and size b with Uq » T 2 /2gb 2 we have 
(see, e.g. [11]) £ ~ T / y/2gU and t\ ~ j^e~ 2lTa ^. Upon using only the first two terms to 
evaluate ( |3.15| ), we find 

.... Airati 

(j(v)) = ^ f ±smv, (3.18) 



and it follows from Eq. ( |3.9a|) that the mean square winding number behaves as 

(W 2 ) = 2hL/T 

~e- 2 ^L. (3.19) 

Although the variance is still proportional to L, the coefficient 2t±/T which multiplies L 
vanishes exponentially fast as a — > oo. 



C. A Single Polymer on a Thick Cylinder 

In the previous section we considered the statistics of winding numbers associated with 
the motion of a polymer on a thin cylinder or, equivalently, a random walker on a ring. For 
a more general cylinder of outer radius R, there are two extreme cases, depending on the 
ratio of the cylinder width R — a and the typical transverse distance the polymer crosses 
in "time" L, namely \J LT / g. In the previous section this ratio was zero. Suppose now 
that this ratio to be small, but nonzero, i.e. < Si^t << i. when using eigenfunctions 
of 7i to evaluate winding numbers, this condition implies that the partition function has 
significant contributions only from one radial mode, the one contributing the lowest energy. 
In this limit, the density profile in the radial direction can be approximated by a constant, 
and an analysis similar to that following Eq. (|3.14j ) leads, for a disorder-free sample in the 
limit L — > oo, to 

J T R 

g7r{R — a ) a 
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The opposite limit, in which the cylinder width is essentially infinite, was considered, also 
in the absence of disorder, in Refs. Q and ||10|| . Here, we confine ourselves to a simplified 



discussion of (W 2 ) using the quantum formalism. According to Eq. ( |3.2| ), the variance (W 2 ) 
can be written, 

« i/2 > = r ^M^^ i[r(TM H ■ <3 - 2i) 

Upon approximating the average in the integrand, 

(±M±m Mr(T)]AMT)] ^ „ l±M. d l^ { a,Ht)]amt')}) 

2T 

^—5{r-r'){A 2 [v{r)]), (3.22) 



we have 



9T r L r 

(W 2 ) = — dr d 2 rV{r,T)A 2 {r), (3.23) 
g Jo J 



where V(r, r) is the probability of finding the random walker at position r at height r. For 
an infinitely thick cylinder, we have P(r, r) = -^f-e~ 9T 1 l 2Tr and (W 2 ) oc (logL) 2 . The 
thermal averages in Eq. ( p.9a| , |379b|) now include many eigenstates. 



The effect of disorder can also be understood by a similar analysis. Again, disorder 
introduces a localization length £, characterizing the eigenstates of the Hamiltonian H(v) 
As long as £ is larger than the cylinder's size, the effect of the disorder is weak. When £ is 
smaller than the cylinder's outer radius, but larger than 2ira, its inner circumference, the 
effective outer radius of the cylinder becomes £. Points on the plane whose distance from 
the obstacle is larger than £ support mostly eigenstates that are insensitive to v. Thus, they 
do not contribute to the winding number. Particles that start their way at such points never 
make it to the hole. Therefore, the winding number for such a system becomes, 



f (logL) 2 for ^<e, £>2na, 



(W 2 ) oc 



(3.24) 

[^log(J) for f >e, £>2vra 
Finally, when £ < 2na the contribution to the winding number comes from a 1-D strip 

around the hole, and the 1-D result applies, namely (W 2 ) oc e~ 27Ta ^L. To summarize, 
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for an infinitely wide cylinder in the long-time limit disorder makes the winding number 
increase, for particles that start their random walk close enough to the hole. In the absence 
of disorder the particles diffuse away from the hole, and thus end up with a small rate 
of winding. Disorder prevents them from diffusing away from the hole, and increases the 
winding number. 



IV. WINDING STATISTICS OF MANY POLYMERS 

We now turn to discuss a system with many polymers (or flux lines). As discussed in 
section 2, the two-dimensional Hamiltonian which describes the physics is that of interacting 
bosons, interacting both mutually and with randomly placed pinning centers. The pinning 
centers we discuss here are assumed to result from columnar defects, and thus create a 
r-independent potential. 

The Hamiltonian is a generalization of Eq. fl2.13| ), namely, 



T 



2 N 



■2g 



1. 



2 N 



^=7r-E -V, + ^A(r,) + £tf(r i ) + 5>(|r i -r i |) (4.1) 



j=l i>j 

where V is the interaction potential and U includes both disorder and the interaction with 
a repulsive cylinder centered on the origin. Note that neither V nor U depend on u, which 
determines only the coupling of the bosons to the Aharonov-Bohm flux. The phase diagram 
of the Hamiltonian (|4. 1|) includes several phases, most notably a superfluid phase and an 
insulating bose glass phase in the limit L — > oo. The properties of these phases are reflected 
in the winding number distributions. The winding number for many directed lines is given 
by Eq. (|2.17| ), and the many particle generalization of Eq. ( |3.9b|) is 



TL r r L r L 

a: 

xA a {r)A p (r'). (4.2) 



{W) = ~J d r < n W> i4 W - Jo dT J dT> I dr I d 2 r'{T\j a {v,r)jp{v'y)]) v ^ 



Several comments are in place regarding Eq. (|4.2|) . First, the density operator n(r) = 
^2i6(r — rj) describes the total density of particles (rather than, e.g., the superfluid density). 
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Second, the current-current correlation function should be calculated for v — 0, i.e., the 
current operator is j Q ( r ) = (1/2(7) Z)i[<K r — r i)Pa + Pafi( r — r j))- Finally, the interpretation 
of the second term becomes clearer upon Fourier transformation. After Fourier analysis, 
the imaginary-time current-current correlation function appearing on the right-hand side of 
Eq. ( f4.2[ ) becomes the linear response function Xapisii w = 0) defined by j a (q, oj = 0) = 
J2/3Xap(ci,u = 0) Ap(q,(jj = 0), i.e., the function describing the current response to a weak, 
time-independent vector potential p5 . 



The current-current correlation function Xapij, T ) = (T[j a (r, r)j / 3(0, 0]) in Eq. (|4.2|) may 
be decomposed in longitudinal and transverse parts. The corresponding Fourier decompo- 
sition reads 



1 



Xo/s(r,r) = —J2 



iqr— ilut 



I \1a% . , sir 1a% 

x\\ (q, ^—^r + x±(q> w ) I <W - 



(4.3) 



However, the vector potential A(r) = for winding around an obstacle is purely trans- 
verse, so only x_|_(q, cu) contributes to the winding number. Upon passing to the limit of 



large system dimensions in the x|/-plane, and rewriting the second term of ( |4.2| ) in Fourier 
space, we have 

(W 2 ) = —j d 2 r(n(r))A\r) -Lj j^ 2 xM^ = 0)|A(q)| 2 (4.4) 
where A(q) is the Fourier transform of A(r). 

A(q) = (4.5) 
\Q\ 

Suppose the directed polymers are confined in an annulus of outer radius R » a, where 
a is the radius of the cylindrical obstacle. As R —>■ oo, we can neglect any distortion of 
the line density due to the obstacle and approximate (n(r)) by its average value, n, far 



from the inner boundary. The first term of Eq. ( [4.4| ) then behaves like -^n\n(R/a), i.e., it 
diverges logarithmically with coefficient proportional to the average density of lines. Because 
|A(q)| 2 = 1/q 2 , the integral in the second term is dominated by the behavior of X-l(Qi u> = 0) 
in the limit q — > 0. Upon imposing upper and lower Fourier cutoffs of a -1 and -R -1 , we again 
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find a logarithmically diverging contribution to the winding number, and the mean square 
winding number in the large R limit is 



TT T 

(W*) = —nHR/a)-- 

ll\g Z7T 



\n(R/a). (4.6) 
However, the transverse response function X-l(<2S w ) is wen known to be related to the normal 



density n n of the equivalent system of interacting bosons [16 



T 

lhnx±(q,o; = 0) = -n n . (4.7) 

q^O g 

Upon defining the superfluid number density n s = n — n n , we are led to our final result, 
namely 

TT 

(W 2 ) = —n s HR/a). (4.8) 
2irg 

Equation (|4.8|) is an exact relation, valid in the limit R » a, between winding number 
fluctuations and the superfluid density of the equivalent boson system. When reexpressed 
in terms of the boson parameters of section 2 it reads 

which is similar in some respects to the Polio ck-Ceperley result (|2.9|) for boson world lines 



winding around a torus ||26|| . Unlike the Pollock-Ceperley formula, however, Eq. ( |4.8|) 
applies directly to a real physical system, namely "directed polymer melts" composed of 
vortex lines in Type II superconductors with columnar defects. As discussed in section 2, 
(W 2 ) determines the linear response of the net winding number to a current through the 
obstacle. 

According to Eq. ([4.8|) the winding number fluctuation around an obstacle for many 
interacting directed polymers is predicted to diverge linearly with the length of the obstacle 
and logarithmically with the cross-sectional area, as one might guess by multiplying the 
result ( |3.20| ) for a single polymer in a moderately thick cylinder by N. The meaning of the 
coefficient becomes clearer if we first define a "thermal deBroglie wavelength" At for the 
polymers by 
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i.e., the directed polymer analogue of Eq. (|2.7f ) with the usual identifications h — > T, 



h/T — > L and m — > g. This length is a measure of the transverse wandering distance of the 
polymers as they traverse the sample. Equation (f4.8|) then becomes 



(r) = (|) 2 n s ln(E/a), (4.11) 

showing that only a part of the "superfluid fraction" n s of lines, i.e., those within a transverse 
wandering distance of the obstacle, contribute to the winding number fluctuations. Corre- 
lated disorder along the f axis should decrease the winding number fluctuations. Equation 
(EHUD shows that this reduction is given entirely by the corresponding reduction in the super- 
fluid density. The reduction in n s for directed lines subjected to various kinds of correlated 



disorder is calculated explicitly in Ref. [23 . 
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FIGURES 

FIG. 1. Slab of Type II superconductor with a dense array of columnar pins confined 
to a cylinder with a dilute concentration of columnar defects outside. Vortices inside the 
cylinder are in the Bose glass phase, while those outside are free to move and entangle in a 
flux liquid. 

FIG. 2. Effect of a current through a cylindrical obstacle on a flux liquid. The lines 
wind around the obstacle with a preferred chirality when the current is on. 

FIG. 3. Boson trajectories in 2 + 1 dimensions projected onto the xy-p\ane. Three 
trajectories with characteristic interparticle distance a and four images of the periodic box 
of size D are shown. The periodic boundary conditions in the time-like direction mean 
that only small ring trajectories appear in the high temperature image at left. Trajectories 
which cross between the periodic cells and lead to nonzero winding numbers are shown in 
the low-temperature image on the right. 

FIG. 4. Mean square winding number (W 2 ) for bosons in 2 + 1 dimensions a function 
of (5 — 1/ksT. A Kosterlitz-Thouless transition occurs for (5 = f3 c . Dashed line is the 
corresponding quantity for directed polymer melts. 

FIG. 5. Trajectories for directed polymer melts with "open" boundary conditions at the 
top and bottom oif a slab of thickness L. 
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